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OPTIMUM  FILTERING  REQUIRED  FOR  BROADBAND 
DETECTION  OF  RANDOM  PROCESSES 

INTRODUCTION 

It  has  long  been  believed  that  optimum  filtering  for  broadband  detection 
consists  of  using  an  EcKart  filter.  The  Eckart  filter  optimizes  the 
deflection  criteria;  however,  it  does  not  optimize  the  criteria  of  usual 
interest,  which  is  maximum  probability  of  detection  (PD)  for  a  fixed  false 
alarm  rate. 

In  this  report,  we  present  the  derivation  of  the  optimum  filter  that 
maximizes  Pq  for  a  fixed  false  alarm  rate.  We  show  by  example  that  the 
optimum  filter  out  performs  the  Eckart  filter  and  other  filters. 

Additionally,  we  present  a  method  of  estimating  performance  when  the  input 
spectra  are  unknown  and  must  be  estimated. 

PROBLEM  DEFINITION 

The  configuration  of  interest  is  depicted  in  figure  1.  The  input  x(t)  is 
composed  of  stationary,  zero-mean,  Gaussian  signal,  s(t),  and  noise,  n(t),  or 
noise  n(t)  alone.  The  sampler  and  summer  effectively  accumulates  M 
statistically  independent  samples  during  an  observation  time  of  T  seconds. 


Figure  1.  Processing  Configuration 

If  the  random  variable  z  exceeds  the  threshold  A  when  signal  is  present, 
a  detection  occurs.  If  z  exceeds  the  threshold  when  signal  is  absent,  a  false 
alarm  occurs.  The  problem  is  to  find  the  optimum  filter,  HQ(f),  that 
maximizes  PD,  for  a  given  probability  of  false  alarm,  Pp,  and  a  given 
observation  time,  T. 
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PROBLEM  SOLUTION 

We  first  define  the  signal-to-noise  ratio  (S/N)  at  the  input  of  the 
square  law  detector  in  figure  1  that  is  required  to  detect  a  signal  with  given 
PD  and  Pp,  i .e. , 


S/N* 


§  X(PD)  -§  1  (Pp) 


(1) 


M  s  BT, 


(2) 


where 


x 

$(X)  -  J  dt ( 2-rr) -1/f 2  exp(-t2/2)  .  (3) 

—  o© 

The  derivation  for  the  above  is  given  by  reference  1.  The  filter  output  S/N 
and  bandwidth  (B)  are  given  by 


and 


ob 


1 

2  Gs(f) 

fdf  1 

[h  (f ) 

|2yf> 

5 


(4) 


(5) 


Equation  (1)  was  derived  by  approximating  the  random  variable  z  as  Gaussian; 
this  is  a  good  approximation  for  M  »  1. 


Substituting  (2)  into  (1)  and  rearranging,  we  obtain 


r 


$ 


S/N  YbT  +  |-1 
1  +  S/N 


(6) 
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Since  <£,  as  given  by  (3),  is  a  monotonical ly  increasing  function,  PQ  can  be 
maximized  by  maximizing  its  argument 


S/NVBT,+  1'1  (Pr) 

- T-TS7N - -  •  (7) 

For  the  filter  output  S/N  <<  1,  as  is  usually  the  case,  and  recognizing  that 
( Pp )  and  T  are  constants,  PD  can  be  maximized  by  maximizing  Q,  where 

Q  —  (S/N)2  B  .  (8) 


Appendix  A  provides  the  derivation  for  the  optimum  filter  |HQ(f )j  2 
that  maximizes  Q.  } HQ ( f )j  2  is  given  by  the  following  equation: 


Vf>l2  - 


"(f)  sN(f) 


31,  -  (91,  -  8I.I-) 

R(f)  — - - 2-2- 


1/2 


(9) 


where 


WQ  is  the  largest  region  where 


R(f)  =  G$(f)/GN(f), 
InF  fdfRn-2(f). 

M 


91 


1  ->  ^o1?' 


Substituting  (8)  into  (6)  and  utilizing  S/N  <<  1,  we  have 


(10) 
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COMPARISON  WITH  OTHER  FILTERS 

Other  filters  are  presently  in  use  for  broadband  processing.  These 
include  the  ideal  passband  filter,  the  Eckart  filter,  and  the  prewhitening 
filter.  These  filters  are  popular  because  of  their  ease  of  implementation  or 
particular  performance  measure  that  they  optimize. 

We  now  compare  the  performance  of  the  optimum  filter,  HQ(f),  with  the 
performance  of  the  forenamed  filters.  The  best  way  to  make  this  comparison  is 
by  a  numerical  example.  Consider  the  following: 


Gs(f)  = 

r.01  exp  (-f/750) 

f  < 

5000^; 

l  o 

f  > 

5000 

S(f)  - 

exp (-f/ 1000) 

for 

O 

A 

M- 

with  integration  time 


T  =  250  sec, and  PF  =  10  6.  (12) 

Equations  (8)  and  (10)  are  used  to  determine  performance.  The  results 
are  presented  in  table  1;  it  shows  that,  for  this  example,  HQ(f)  yields  the 
best  PD. 

Table  1.  Filter  Comparison 


F i 1  ter 

Power 

Transfer 

Function  [h ( f  )|  p 

Q 

PD 

Optimum 

Equation  (9) 

0.112199 

0.71 

Eckart 

GS(0/G2  (f) 

0.098704 

0.58 

Prewhitening 

0.096571 

0.56 

Ideal  Passband 

1 

0.085089 

0.44 
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If  filtering  is  performed  in  the  time  domain,  HQ(f)  is  no  more 
difficult  to  implement  than  the  other  filters.  However,  it  is  more  difficult 
to  design.  The  same  holds  true  if  filtering  is  performed  in  the  frequency 
domain,  except  that  the  ideal  passband  filter  is  trivial  to  implement.  The 
other  advantage  to  the  optimum  filter  is  that  the  whole  signal  band  need  not 
be  processed.  The  cutoff  for  the  optimum  filter  in  this  example  is  3455  Hz, 
whereas  the  other  filters  cut  off  at  5000  Hz. 

UNKNOWN  SPECTRA 

When  G«.(f)  and/or  G^(f)  are  not  completely  known,  then  PD  is  a 
function  of  the  random  variable  Q.  Two  characteristics  of  PQ  that  are  of 
interest  are  its  expected  value  and  variance. 

Writing  (10)  as 


PD=F(Q),  (13) 

and  making  use  of  section  7.7  of  reference  2,  we  obtain  the  following 
approximations : 

E[Pd] ~  F ( E[Q] )  +  F"(E[Q])  Var[Q]/2, 

Var[PD]  *  [F* (E[Q])]2  Var[Q]  .  (14) 

Substituting  (3)  and  (10)  into  (13),  there  follows 
F ' ( E [Q ] )  =  (2 u)-1/2  T/2(E[Q]T)-1/2  exp[-[(E[Q]T)1/2  + §_1(PF)]2/2], 

F"(E[Q])  =  —  ( 2ir )  1  / 2  T2/4(E[Q]T)-1  [( (E[Q]T)1/2  +$-1(Pp))  +  (E[Q]T)-1/2] 

*exp[-[(E[Q]T)1/2  +§-1(Pf)]2/2].  (15) 

Appendix  B  is  used  to  evaluate  E[Q]  and  Var[Q]  (see  (B-8)). 
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Example  , 

Consider  the  example  from  Appendix  B  with 

T  =  450  sec,  Pp  =  io"7,  5_1(PF)  =  -5.2  .  (16) 

Since  Q(R)  =  0.08,  using  (10),  we  obtain 

PD(R)  =  0.7881  .  (17) 

Since  E[Q]  =  0.08  and  Var[Q]  =  0.000256,  substituting  into  (14)  and  (15),  we 
obtain 


F'(E[Q])  =  10.863433, 

F"(E[Q])  =  -393.799455,  (18) 


and  therefore 


E[PQ]  =  0.7377, 

Var  [PD]  =  0.0302.  (19) 

The  error  in  the  estimate  of  the  signal  spectrum  has  caused  a  decrease  in 
expected  performance. 


CONCLUSION 


The  derivation  for  the 
detection  for  a  given  false 
improvement  can  be  achieved 
be  degraded  when  the  signal 


broadband  filter  that  maximizes  probability  of 
alarm  rate  has  been  provided.  Performance 
using  this  filter.  The  detection  performance  will 
or  noise  spectrum  must  be  estimated. 
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APPENDIX  A 

DERIVATION  OF  OPTIMUM  FILTER  POWER  TRANSFER  FUNCTION 

For  notational  convenience,  we  drop  the  f-dependence  in  (4)  and  (5) 
temporarily.  Then  defining 


A  ,  |h|2Gn,  r  =  gs/gn, 

we  have  to  maximize,  by  choice  of  A,  the  quantity 


Actually,  we  will  consider  the  slightly  more  general  form 


( A— 1 ) 


( A— 2 ) 


(A-3) 


where  R^  and  R2  are  given  known  non-negative  functions  of  frequency,  f. 

The  integrals  in  (A-3)  are  over  the  band  W  of  interest,  and  would  be  written 
more  accurately,  for  example,  as 

a  =  JaRj  =  jdf  A(f )Rj(f )  .  (A-4) 

W 


Notice  that  the  function  A  defined  in  (A-l)  can  never  be  negative. 


To  determine  the  optimum  A,  we  let 


A(f)  =  A0(f)  +  en(f)  ,  (A-5) 

where  AQ  is  the  optimum,  and  n  is  an  arbitrary  perturbation.  Substituting 
(A-5)  into  (A-3),  we  obtain 


A-l 
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Q0  +  «Q- 


[f(A0  +  er,)R^] 

r 

[h  * ««)' 

2 

f(A  o  +  en)2R2 

(A-6) 


Taking  a  partial  of  (A-6)  with  respect  to  e,  and  then  setting  e  =  0,  we  get 

[So  Yo  4  ao  j*n  R1  -  ao  (2Bo  yo  fn  +  Bo  2  J\  n  R2}]/(fo  Yo)  =  0  for  any  n*  (A-7) 

Cancelling  irrelevant  constants,  none  of  which  can  be  zero,  we  find  that 


Yo  1 

Rj  -  j-  -  R2  Aq j  =  0  for  any  n.  (A— 8) 

The  only  way  this  can  be  is  if  the  bracketed  term  is  zero  for  all  f  e  W.  Thus 
re-substituting  the  f  dependence,  the  unique  optimum  A-function  satisfies 


R2(f)  Ao(f)  = 


for  f  e  W 


(A-9) 


Now  if  R2<f)  were  zero  in  some  subregion  of  W,  where  R ^ (f )  is  not  zero,  it 
is  obvious  from  (A-3)  that  Q  can  be  made  equal  to+«>by  simply  choosing  A(f)  to 
be  non-zero  only  in  that  subregion,  for  then  y=0.  Therefore,  we  exclude  this 
case  and  require  that  R2(f)  >  0  for  f  e  W. 


Also,  since  A  (f )  can  never  be  negative,  we  must  replace  any  negative 
values  that  arise  in  (A-9)  by  zero.  Accordingly  we  define  region  WQ  such 
that 


Rl(f)>2T  for  f  '  Wo 


Then  the  optimum  A  function  is 


A0(f)  - 


2y 

O  1  p  /  f  \  o 

^r~^rrj  Ri(f)  -w: 


for  f  e  W„ 


A- 2 


0  for  remaining  f  e  W 


(A-10) 


( A— 1 1 ) 
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Several  useful  observations  should  be  made  at  this  point.  First,  from 
( A-3 ) ,  the  absolute  scale  of  A  or  AQ  does  not  affect  Q;  so  we  will 
eventually  be  able  to  discard  the  scale  factor  2yQ/a0  in  (A-ll).  However, 
we  must  temporarily  retain  it  in  order  to  solve  for  the  constants  ctQ,  s0, 
yq  defined  in  (A-3). 

Second,  we  do  not  need  yQ  at  all  since  it  is  merely  a  scale  factor  in 
(A-ll);  and  we  do  not  need  to  know  aQ  and  eQ  separately,  but  only  their 
ratio.  Also,  aQ,  e0,  y0  must  all  be  positive  since  they  are  integrals 
of  positive  functions. 


Finally,  the  form  of  the  optimum  A  ff)  in  (A-ll)  indicates  that  if 

AQ(f )  is  zero,  it  is  zero  in  the  frequency  regions  where  Rj(f)  is  minimum, 

regardless  of  what  R£(f)  is.  This  is  a  very  useful  observation  for 

determining  W  . 

3  o 

Now  we  substitute  (A-ll)  in  the  individual  terms  of  (A-3)  to  determine 
the  constants  aQ,  e0,  Yq.  Dropping  the  zero  subscripts  on  a,  e,  y  for 
notational  convenience,  we  have 


where 


I 


n  ' 


f  R?(f) 

]  df  ^rry 

Wo 


( A— 13) 


A-3 


Simplifying  (A-12),  we  have  the  three  simultaneous  nonlinear  equations 


(A-14) 


At  first  sight,  solution  of  this  set  looks  formidable  indeed.  However, 
observe  that  the  third  equation  is  not  independent  of  the  first  two;  in  fact, 
2b  times  the  first  equation  minus  a  times  the  second  equation  yields  the  third 
equation.  So  we  discard  the  third  equation  in  (A-14). 

Next,  eliminating  the  unwanted  y  from  the  first  two  equations  in  (A-14), 
we  obtain  (restoring  the  zero  subscripts  on  a,  b,  y) 


9 


( A— 15 ) 


which  is  an  equation  in  exactly  the  one  unknown  we  fundamentally  need  in 
(A-ll).  The  solutions  are 


(A-16) 


Define  square  root 


1/2 


(A-17) 


In  order  that  a  real  solution  for  ctQ/B0  exist,  as  required  by  definitions 
(A-3),  we  must  also  have  discriminant 
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Thus  when  we  pick  the  region  WQ  in  (A— 13) ,  it  must  be  done  with  constraint 
(A-18)  in  mind.  If  (A-18)  is  not  satisfied,  the  corresponding  choice  of  WQ 
is  disallowed. 


It  should  be  observed  that  (A-18)  is  a  very  tight  restriction  on  Iq, 
Ip  12*  For  when  we  couple  it  with  the  Schwartz  inequality  (see  (A-13)) 


we  have 


8  T  T  Z 

9  0  2  -  1  -  0  2 

In  order  to  determine  which  of  the  two  positive  roots  for  aQ/eo  is 
the  correct  one  to  retain  in  (A-16),  we  must  substitute  it  into  Qq  and  find 
which  yields  a  maximum.  From  (A-16)  and  (A-17), 


s) '  (r*v)  *  th  (18  " 


8  IQ  I2  ±  6  I 


is)  • 


(A— 19) 


and  from  the  first  equation  of  (A-14), 


T0-  2  '2  ’2  -1 


4  i0  i2  -  3  ir  + 


X*  ‘Is) 


( A-20) 


Then  from  (A-3),  (A-19),  and  (A-20), 


2  2 

0  -  -2.  _o  _ 

g°  "  q2  rn  ” 

e0  To 


2  i; 


36  IQ  if  I2  “  8  lQ  lz  ~  27  1 


i  +  'i  s3) 


(A— 21 ) 


This  is  obviously  maximized  by  choosing  the  lower  sign  here  and  in  (A-16).  So 
we  have 


A-5 
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3  Ix  -  S 

2  17 


( A— 22 ) 


and 


'°-7?A 


36  IQ  if  I2  -  8  Iq  l\  -  27  ij  +  I 


.>•) 


2  i; 


27  I 


1  (!0  '2  -  ‘0 


+  IQ  I2  D  +  Ix  D  S 


( A— 2  3 ) 


All  the  terms  in  (A-23)  are  obviously  positive,  the  Ig  term  following 

from  Schwartz's  inequality  under  (A-18). 

We  can  now  drop  irrelevant  scale  factors  in  (A-ll)  and  state  that 


Vf>  =  T^ITT 


l 


Rl(f)  ~  7F" 


Rj_ (f )  - 


3  ll  -  S 


4  I 


0  J 


for  f  e  W„ 


(A-24) 


The  quantities  Ig,  1^,  l2,  and  S  needed  in  (A-22)-(A-24)  are  given  by 
(A-13)  and  (A-17).  We  must  keep  in  mind  that  (A-24)  can  never  be  negative  and 
that  (A-18)  must  be  satisfied.  The  optimum  filter  power  transfer  function 
follows  from  (A-l)  as 


If 


frequencies 


were  known,  we  could  easily  determine  WQ  by  finding  the 

a  a 

where  R^(f)  =  >  as  noted  earlier  under  (A-ll);  but 


( A— 25 ) 


itself 


A-6 
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depends  on  WQ  through  the  quantities  IQ,  1^,  I^.  The  way  to  determine  y—  is 

OL 

as  follows.  Hypothesize  a  value  for  ;  call  it  c.  Determine  Wc  from  R^(f) 

by  solving  ( f c )  =  c,  and  then  choosing  Wc  as  that  region  where 

R^f)  >  c;  there  may  be  several  values  of  fc.  (There  has  been  no  need  to 

consider  R^( f )  up  to  this  point.)  Compute 

r  Rj(f) 

‘n  '  J  df  1CTTT  f0r  "  "  '>  2  *  (A-26) 

w  d 

c 

which  depend  on  c.  Compute  D  from  ( A— 18) ;  if  D  <  0,  the  hypothesized  c  value 
is  disallowed.  Otherwise,  compute  S  =  Yd1,  and  then  a/8  from  (A-22).  Now  we 
will  have  the  correct  value  of  c  when  we  get  back  a/ (2b)  =  c;  so  we  have  the 
identity 


c 


h  -S 

IT 


3  I 


I"  (9 


"  8 


tt: 


■r 


(A-27) 


or 


2  IQ  c2  -  3  Ij_  c  +  I2  =  0  .  (A-28) 

This  last  equation  is  identical  to  (A-15).  Solution  of  (A-27)  or  (A-28) 
(numerically  or  analytically)  determines  <*0/(2b0).  Negative  values  for 
AQ(f )  are  never  encountered  or  contemplated  via  this  approach.  The  only 
check  to  make  is  that  D  is  not  negative;  if  it  is,  the  corresponding  selection 
of  c  is  disallowed  and  must  be  modified. 
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Example  1 


We  consider  first  the  special  case 


R^f)  =  R(f),  R2(f)  =  R2(f), 


where  R(f)  is  a  two-level  function,  i.e.. 


R(f) 


for  0  <  f  <  1 


for  1  <  f  <  2 


( A— 29 ) 


( A-30 ) 


This  example  has  equal  bandwidths  for  each  of  the  two  frequency  intervals, 
namely  (0,  1)  and  (1,  2),  normalized  to  1  Hz  for  convenience.  Without  loss  of 
generality,  we  let  r1  <  r2  and  define 


r 


r  >  1 


(A— 31 ) 


We  now  employ  the  procedure  described  in  (A-26)-(A-28) .  First  assume 


c  <  r 


1 


(A-32) 


Then  Wc  =  (0,  2),  and  from  (A-13),  (A-29),  (A-30),  and  (A-18), 

In  =  £  df  Rn_2(f )  =  rj”2  +  r£-2  =  q2~"  +  q2‘n  ,  (A-33) 

0 

where  qk  =  l/rk.  Then 

1/2 

!0  =  qi  +  <*2*  !1  =  ql  +  q2’  12  =  2’  S  =  (18  ql  q2  "  7  ql  -  7  q2)  •  (A“34) 

The  identity  (A-27)  yields  immediately 
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c  = 


/  9 

3  q,  '  3  q2  -  (18  qt  q2  -  7  qt  -  7  q2)  3  +  3  q2  -  S 


W® 


Wi) 


(A-35) 


for  the  value  of  jg-,  provided  only  that  the  discriminant  is  non-negative. 
Satisfaction  of  this  requirement  yields  (in  conjunction  with  (A-31)) 


Q  +  (]aF? 

1  £  r  <  -  j  -  ■<c  =  2.0938 


( A— 36 ) 


Also  we  have  to  confirm  that  (A-32)  is  true;  i.  e.. 


3  qx  +  3  q2  -  S  ? 


W ?) 


iri  =  q^ 


( A— 37 ) 


This  can  be  manipulated  into  the  form 


4  q2  +  (2  q2  -  ql)  *  '’l  S  i  0  ' 


(A-38) 


But  since  this  is  obviously  always  true,  the  only  restriction,  for  this  case 
of  assumption  (A-32),  is  (A-36).  That  is,  (A-32)  is  an  allowable  assumption 
provided  that  (A-36)  is  satisfied. 

We  now  consider  the  complementary  case  to  (A-32);  namely,  we  assume 


c 


(A-39) 


Then,  Wc  =  (1,  2) ,  and 


2 

>n-Jdf  Rn-2(f) 
1 


r 


n-2 

2 


(A-40) 


A-9 


TR  No.  6999 
that  is. 


2  2 
Iq  =  q^,  =  q2>  1 2  =  1 >  D  =  q2  >  0,  S  =  q2 

Identity  (A— 27 )  yields  immediately 


for  the  value  of  aQ/(2B0).  In  order  to  satisfy  (A-39),  we  must  have 


(A-41) 


(A-42) 


2 

2~  >  rl  or  r  >  2  .  (A— 43) 

We  now  observe  that  conditions  (A-36)  and  ( A— 43 ).  cover  all  the  possible 
values  of  r  in  (l,+«>);  in  fact,  they  overlap  in  the  interval  (2,  2.0938).  In 
order  to  determine  whether  to  select  (A-35)  or  (A-42)  in  this  overlap  interval 
(and  both  are  allowed),  we  must  evaluate  Qq  in  this  interval  and  take  the 
larger  value. 


We  find  that,  from  (A-23)  and  (A-34), 


^o  =  *2 


2 

13(1  +  r^)  +  36r(l  +  r^  -  3r)  +  (1  +  r)  (l8r  -7-7  r^) 


3/2 


2  (i  *  r2)' 

for  condition  (A-36);  while  from  (A-23)  and  (A-41), 


%  -  4 


for  condition  (A-43).  We  also  find  that  Qq/ r ^  in  (A-44)  decreases 

monotonical ly  from  the  value  2  at  r  =  1,  to  the  value  1  at 

r  =  2.06353  . 


(A-44) 


( A— 45 ) 


(A-46) 
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Thus,  for  r  <  r  ,  the  maximum  Qq  is  given  by  (A-44),  while  for  r  >  r  , 
the  maximum  Qq  is  given  by  (A-45).  Q0  is  continuous  with  r,  but  has  a 
discontinuous  slope  at  r  . 


As  for  the  optimum  AQ  function,  we  find,  by  use  of  (A-24)  and  (A-35), 


that  when  1  <_  r  £  2.0938,  with  q  =  (18^^  -  7r2  -  7r2)^2, 


Vf>  - 


for 


(A-47 ) 


Since  the  absolute  scale  of  AQ(f)  is  not  important,  only  the  ratio  of  filter 
gains  is  essential.  It  is,  after  considerable  manipulation. 


(?)  1/2 
Ao  ^  1  -  r  -i-  (  18r  -  7  -  7r2) 

A^(f)  "  -  rj 


( A-48 ) 


(When  r  =  2,  this  ratio  is  3/2  via  a  limiting  operation.)  For  r  =  1,  the 
ratio  in  (A-48)  is  1  as  expected,  since  the  function  R(f)  in  (A-30)  is  flat 
over  (0,  2).  At  the  other  extreme,  where  r  is  the  largest  allowable  value  for 
this  case  (see  (A-36)),  the  ratio  of  filter  gains  is 


Ao2)(f>  11  +  6  V? 

iPm  =  __7 — 


2.7836 


for 


r 


9  + 


4V7 

7~ 


(A-49) 


while  at  r  =  r  =  2.06353,  the  ratio  is  1.8440.  Thus,  ratio  (A-48)  increases 
monotonical ly  from  1  (at  r  =  1)  to  the  value  in  (A-49).  Observe  that  this 

A^(f)  is  not 


ratio  is  not  +«t>  at  this  upper  limit  for  r  in  this  case;  i.e., 
zero  at  this  value  of  r  =  (9  +  4YT)/7. 
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For  the  complementary  case,  (A— 43) ,  we  use  ( A— 24 )  and  (A-42)  to  get 


Vf>  - 


for 

f  e  (o,  1) 

for 

f 

f  e  (1,  2) 

J 

A 

1 

Iv 

-s 

,  while  (A-47) 

(A— 50) 


imum 


for  1  <  r  <  r  So  at  r  =  r  ,  A^(f)  jumps  abruptly  from  a  positive  value  to 


zero. 


Example  2 


Here  we  consider  continuous  functions  of  frequency 


Rj_(f)  =  exp(-Tf)  for 


f  >  0  I 


f^2  ( f )  =  exp(-axf)  for  f  >  0 


W  =  (0,+«*)  ,  .  (A-51) 


where  a  is  a  dimensionless  known  constant;  if  a  =  2,  then  R2(f)  = 
R^(f),  the  original  case  of  interest  in  (A-2). 

We  use  the  procedure  described  in  (A-26)-(A-28) .  Let  c=l. 
Then  Wc  =  (0,  fc),  where  fc  is  the  solution  of 


Rl(fc)  =  exP  (-Tfc)  =  c;  rfc  = 

Then, 


I 

n 


f  R?(f)  oc 

J  df  rJJT  =  o 


R21fT 


df  exp  (-ntf  +  arf)  = 


exp[(a  -  n)Tfc]  -  1 
(a  -  n)r 


(A-52) 


( A— 53 ) 
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The  only  case  we  will  consider  here  is  a  =  2;  then,  using  (A-52), 


t  1  1  -  c  T  1  1  -  c  T  £  1 

Iq  ~  T  2c2  ’  !l  "  T  c  ’  *2  "  fc  _  “  r^n(c) 


( A— 54 ) 


Identity  (A-28)  yields 


1  -  c2  -  3  (1  -  c)  -  Jn(c)  =  0  , 

(c  -  1)  (c  -  2)  +  in(c)  =  0  .  (A-55) 

There  are  only  two  solutions  of  this  equation;  they  are  c  =  1  (disallowed)  and 


c  =  .316197  =  .  (A-56) 

po 

Then  (A-52)  yields 

Tfc  =  1.15139  ,  .  (A-57) 

The  value  of  Qq  is  obtained  by  substituting  (A-54)  and  (A-56)  in 
(A-23),  yielding 


TQ0  =  .373997  .  (A-58) 

The  optimum  filter  follows  upon  use  of  ( A— 51 ) .  (A-54),  and  (A-56)  in  (A-24); 
i  .e. , 


rU  -  r 

A0(f)  =  — 2-u—  for  0  <  u  <  1  ( A-59) 

c 

where 

u  =  f-  .  ( A-60 ) 

c 

A  plot  of  (A-59)  is  given  in  figure  A-l.  The  optimum  AQ(f)  goes  to  zero  at 
f  =  fc;  it  also  lends  heaviest  emphasis  to  frequencies  near  .4fc,  rather 
than  at  zero  where  R^(f)  and  R£(f)  are  largest. 
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<4 


Figure  A-l.  Optimum  AQ(f )  of  (A-59) 
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APPENDIX  8 


SENSITIVITY  ANALYSIS  OF  FILTER 


It  often  happens  that  Gs  and/or  GN  are  not  completely  known  but  must 
be  "guessed"  or  estimated.  Furthermore  H(f)  usually  cannot  be  realized 
exactly.  All  of  these  factors  influence  performance  and  a  derived  optimum 
filter  may  not  be  optimum  in  fact. 

When  R  (the  estimate  of  R)  differs  from  the  true  R  in  a  deterministic 
manner,  Q  can  be  computed  via  (B-l).  However,  if  R  differs  from  R  in  a 
probabilistic  manner,  then  Q  is  a  random  variable  and  one  can  only  compute  the 
expected  performance. 

From  (A-l)  and  (A-3)  in  appendix  A, 


(B-l) 


where  R  =  Gs/gn  ,  A  =  JhJ2Gn 


The  subsequent  analysis  assumes  that  only  G^  must  be  estimated  (i.e., 

Gn  is  exactly  known  and  H  can  be  exactly  implemented).  The  analysis  can  be 
extended  to  the  case  of  implementation  errors  in  the  filter  and  unknown  G^. 

We  assume  that  G^  is  estimated  by  breaking  frequency  region  W  into  n 
equi length  subregions  W^,  1  <  i  £  n,  where  the  power  in  each  subregion  is 
measured.  Thus,  n  must  be  chosen  large  enough  so  that  Gs  and  GN  remain 
fairly  constant  within  each  subregion.  The  problem  is  then  transformed  into 
the  discrete  domain.  Namely,  we  can  now  rewrite  (B-l)  as 


n  I4 


o 


4 


(B-2) 
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where  Q  is  a  function  of  the  random  variables  ,  1  <  i  <  n. 

Q  can  be  expressed  as  a  multidimensional  Taylor  series,  i.e.. 


n 

i=l 


Q(R)  =  Q(B)  +  lEaQ/aRj  (R^) 

I  K 


+  i_ 


1  n  n  9  1 

a^Q/aR.aR.  (r._r.)(r._r  .) 

21  uijfi1  Jk  1  1  j  y 


(B— 3) 


If  (R^-R^)  is  small,  then  E[Q]  and  Var[Q]  may  be  approximated  by 
neglecting  high  order  terms  in  the  Taylor  series.  Furthermore,  if 


ECR-j-Ri]  =  0, 

Var[R.]  =  CT?, 

and  all  the  ^R^  are  independent,  then 


E[Q]  s  Q(R)  +  i  Ji  32Q/aRi2L  » 

d  1-1  1  IR  1 

Var[Q]^  ^  (3Q/9R.)2L  o2.  . 

1-1  1  IR  1 

The  chain  rule  can  now  be  used  to  evaluate  aQ/aR^  and  a^Q/aR?. 
Making  use  of  (B-2),  we  find  that 


a  Q  /  aR,- 1^  = 

|R  B 


4Aia" 


2A?  Ri  a4 


( B— 4) 


(B-5) 
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2  2 
3  Q/aR^ 


R  f 


'3A -  a 


4Ai 


2A] 


2A.2  R.  a3 


4A2  R2  a4 


4A.Ria3  +  a4'' 


Rearranging  (B-6),  we  obtain 


aQ/aRiJ^  =  Q(R)j4Ai/«  -  2A.2  R./y]_  , 

-  =  Q(^)  12A  2/a2  -  16A  3  R./(ay) 
R  L  1  11 

+  8An.4  Rt2/Y2  -  2Ai2  . 


a2Q/aRi2 


Substituting  (B-7)  into  (B-5),  we  obtain 


E[Q]  =  Q(R) 


Var[Q]  =  Q^(R) 


B)jjL  +  \  (12A.2/a2  -  16A.3  R./(ay) 

+  8A.4  R^/y2  -  2Ai2  /y)  o.2^  , 


’2  *'M%  °' 


i=l 


(4A./o  -  2A 


Example 


(B-6) 


(B-7) 


( B— 8) 


Reconsider  Ex.  1  from  Appendix  A  with 
E  [Gsl]  =  2  . 

Var  [G$1]  =  0.04  , 

E  [GS2]  =  2  , 

Var  [G-?]  =  0.12  , 

S  -  10. 


( B— 9 ) 
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Then, 


( B— 10 ) 


For  this  trivial  example,  by  inspection  let 


( B— 1 1 ) 


Then,  from  (B-2) 


A:  =  10 
a2  =  10 


( B— 12 ) 


and  at  R 


8  =  20  , 

Y  =  8  , 

Q  =  0.08. 

Substituting  ( B— 10 ) ,  (B-12),  and  (B-13)  into  (B-8) 

E[Q]  2  0.08  , 

Var  [Q]  sr  0.000256, 

Std  dev  [Q]  s  0.016  . 

So  Q  lies  most  often  in  the  region  0.08  ±  0.016. 
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